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We study the stability of a spherically symmetric perturbation around the flat Friedmann- 
Lemaitre-Robertson-Walker spacetime in the ghost-free bigravity theory, retaining nonlinearities 
of the helicity-0 mode of the massive graviton. It has been known that, when the graviton mass is 
smaller than the Hubble parameter, homogeneous and isotropic spacetimes suffer from the Higuchi- 
type ghost or the gradient instability against the linear perturbation in the bigravity. Hence, the 
bigravity theory has no healthy massless limit for cosmological solutions at linear level. In this 
paper we show that the instabilities can be resolved by taking into account nonlinear effects of the 
scalar graviton mode for an appropriate parameter space of coupling constants. The growth history 
in the bigravity can be restored to the result in general relativity in the early stage of the Universe, 
in which the Stiickelberg fields are nonlinear and there is neither ghost nor gradient instability. 
Therefore, the bigravity theory has the healthy massless limit, and cosmology based on it is viable 
even when the graviton mass is smaller than the Hubble parameter. 

PACS numbers: 04.50.Kd, 98.80.-k 


I. INTRODUCTION 


General relativity (GR) is now widely accepted as the 
low-energy effective theory of gravity and has passed a 
number of observational and experimental tests [l|. It 
can also be viewed as the unique theory of a massless 
spin-2 particle, namely graviton. Although theories with 
a massive spin-2 field are one of the most natural ex¬ 
tensions of GR, such theories have suffered from many 
problems. Fierz and Pauli proposed a massive spin-2 
field theory, which is known as the unique theory of lin¬ 
earized massive gravity free from a ghosty sixth degree of 
freedom Q. The Fierz-Pauli (FP) theory has appropri¬ 
ate five degrees of freedom, and there is no ghost mode 
on the Minkowski background. However, the FP the¬ 
ory has no Newtonian limit even in a massless limit. 
When we consider a naive massless limit in the FP the¬ 
ory, the gravitational behaviour is not restored to the one 
in Newtonian gravity due to the existence of the extra 
gravitational degree of freedom, which is called the van- 
Dam-Veltman-Zakharov (vDVZ) discontinuity Q. Vain¬ 
shtein then proposed that the vDVZ discontinuity can 
be evaded by taking into account nonlinear mass terms, 
and the extra mode is screened inside the so-called Vain¬ 
shtein radius, to recover the standard gravitational inter¬ 
action mediated only by the helicity-2 modes Q. Boul- 
ware and Deser pointed out, however, that such nonlinear 
terms ruin the structure of the FP theory and reintroduce 
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the ghost instability associated with the sixth degree of 
freedom Q. Although this nonlinear ghost, often called 
Boulware-Deser (BD) ghost, appears in any simple non¬ 
linear extensions of the FP massive gravity theory, it was 
shown in 2010 by de Rham et al. that the special choice 
of the mass term can eliminate such a ghost state at the 
decoupling limit Q , and later the proof was extended to 
fully nonlinear orders 


However, this theory, often dubbed de Rham- 
Gabadadze-Tolley (dRGT) theory, has still been prob¬ 
lematic in the context of cosmology. It was revealed 
that the original dRGT theory does not admit any 
nontrivial flat or closed Friedmann-Lemaitre-Robertson- 
Walker (FLRW) universe, if the fiducial metric for the 
Stiickelberg fields is Minkowski’s one [Il|. While this is¬ 
sue can be resolved either by open FLRW solutions 121 
or by repla cing the Minkowski fiducial metric with an 
FLRW one [I 3 , it was later shown that all homogeneous 
and isotropic FLRW solutions in the dRGT theory are 
unstable due to either a linear ghost, called Higuchi ghost 
0 , or a new type of nonlinear ghost 0. The ghost of 
the former type was found present already in the FP 
theory when it is constructed on a de Sitter background. 
The scalar graviton mode in the five degrees of freedom in 
the FP theory becomes a ghost when the graviton mass 
is below the so-called Higuchi bound, which is defined by 
a cosmological constant, leading to the conclusion that 
the theory has no healthy massless limit with non-zero 
cosmological constant. 

On the other hand, when we take a curved fiducial ge¬ 
ometry, e.g. the FLRW or an inhomogeneous background, 
it may be natural to promote it to a dynamical one. In 
fact the dRGT massive gravity theory has been general¬ 
ized to such a theory with two dynamical metrics, thereby 
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called a bigravity theory, which is still free from the BD 
ghost 13 . It contains a massless spin-2 field and a mas¬ 
sive spin-2 field, with total seven degrees of freedom in 
the gravity sector [17 1. 

Phenomenologically, the bigravity theory has many in¬ 
teresting features. One of the biggest mysterious prob¬ 
lems in modern cosmology is the origin of the current 
accelerating expansion of the Universe [l^. Based on 
the bigravity theory, self-acceleration in all types of the 
FLRW universe is allowed [ii-ii. Therefore, the present 
acceleration of the Universe can be explained by adding 
a mass to the graviton (see also |lll4l3l IlSl . I26l429l| about 
the cosmology in the dRGT theory). Furthermore, since 
the bigravity theory contains two metrics, a matter field 
in the fiducial metric sector is naturally introduced, and 
it may act as a dark matter component in the physical 
sector [25l.[^. 

However, although the background dynamics in the bi¬ 
gravity cosmology is viable, a problem still exists at the 
perturbation level 31 -[s^- In the late stage of the Uni¬ 
verse, the perturbations around the cosmological back¬ 
grounds are stable. However, the perturbations suffer 
from either the ghost or gradient instability in the sub¬ 
horizon scale, when the Hubble expansion parameter is 
larger than the effective graviton mass. 

Since the bigravity theory contains a massless spin-2 
field, the bigravity theory always contains GR solutions 
exactly as a special case. If the two background met¬ 
rics are proportional, which we call a homothetic solution 
[ 3 ^, the basic equations are reduced to two sets of the 
Einstein equations in GR with a cosmological constant 
originated from the mass term. That is, in homothetic 
spacetimes, the bigravity reproduces the background evo¬ 
lution of the universe identical to the one in GR. The lin¬ 
ear perturbations around a homothetic solution are easily 
decomposed into two eigenstates: the massless and mas¬ 
sive graviton modes. Here, the massive graviton mode is 
given by the FP theory on a GR solution. The perturba¬ 
tions in this mode suffer from the gradient instability in 
the decelerating Friedmann universe, while they have the 
Higuchi-type ghost in the accelerating one [sg. There¬ 
fore, the cosmological instabilities in the bigravity theory 
are similar to those in the FP theory on a cosmological 
background and are related to the fact that the linear 
massive spin-2 field has no healthy massless limit on the 
curved background. 

However, such an instability is quite obscure in the 
physical interpretation, since the natural expectation 
would be that the massive theory should be restored to 
its corresponding massless theory when the energy scale 
of the background spacetime is higher than the mass. 
Hence, the instability should be resolved without either 
a modification of the theory or an extra ingredient, if the 
bigravity theory is a reliable theory in such an energy 
scale. The instability may simply hint the possibility that 
the linear perturbations are no longer valid. Therefore, 


before we conclude the bigravity theory breaks down in 
the early stage of the Universe, it is instructive to take 
into account the nonlinear interactions. 

For this reason, we will consider cosmological pertur¬ 
bations with nonlinear effects in the bigravity theory. 
As we will see in Sec. IHI B[ the instability arises from 
a scalar graviton mode. Therefore, we will focus only on 
this mode and also restrict our analysis to a spherically 
symmetric configuration, for simplicity, in which a grav¬ 
itational degree of freedom is only given by the scalar 
graviton mode. 

In order to find a stable model of the early Universe 
in the massive gravity theory (or the bigravity theory), 
there are several different approaches. One possibility to 
obtain the viable cosmology with a massive graviton is 
to extend the theory [37H4fll | . This approach may be jus¬ 
tified by the reasoning that the bigravity or the dRGT 
theory may necessarily be modified at a high energy scale 
to realize a massive graviton resulting from spontaneous 
symmetry breaking, and to protect the validity of the 
theory in the UV regime. Another possibility is to intro¬ 
duce a doubly coupled matter field. Although a matter 
field interacting with both metrics suffers from a reap¬ 
pearing BD ghost in generic situations, some ghost-free 
matter coupling was discussed at a low-energy scale or in 
a homogeneous Friedmann background universe Hliil- 

El- 

In the present paper, however, we will show that both 
the Higuchi-type ghost and the gradient instability can 
be resolved by the nonlinear effects of the scalar graviton 
mode in the ghost-free nonlinear bigravity theory without 
either a modification of the theory or an additional mat¬ 
ter held. When the Hubble parameter is larger than the 
graviton mass, we hnd the new cosmological solutions 
in the sub-horizon scale for which the two spacetimes 
are still approximately homogeneous and isotropic, and 
these two foliations are related by a nonlinear coordinate 
transformation. That is, the Stiickelberg helds become 
non-linear, while the spacetime perturbations are quite 
small. 

The paper is organized as follows. Introducing the 
ghost-free nonlinear bigravity theory, we show the bi¬ 
gravity theory contains GR solutions as a special case in 
Sec.lE In Sec, mil we derive the quadratic action around 
the GR solution, which can be decomposed into massless 
and massive graviton modes. We discuss that the FP the¬ 
ory has an instability on the cosmological background. In 
Sec. IIVI retaining nonlinearities of the Stiickelberg fields, 
we study a spherically symmetric perturbation around 
the flat homothetic FLRW background without matter 
perturbations, and we discuss the stability in the early 
stage of the Universe. Both Higuchi ghost and gradient 
instability can be resolved by the nonlinear effects of the 
Stiickelberg fields in the early stage of the Universe. In 
Sec. El we introduce the matter perturbation and study 
stability of the matter perturbations both in the early 
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stage and the late stage of the Universe. We discuss a 
possible transition from GR phase to the bigravity phase 
in the period when the graviton mass is comparable to 
the Hubble parameter in Sec. ED We summarize our re¬ 
sults and give some remarks in Sec. IVIIl In Appendix El 
we present the Vainshtein screening in static spherically 
symmetric configurations with a cosmological constant as 
well as twin matters. 


II. NONLINEAR BIGRAVITY THEORY 


A. Hassan-Rosen bigravity model 


In the present paper, we focus on the ghost-free bi¬ 
gravity theory proposed by Hassan and Rosen [16| , whose 
action is given by 


•S' = ^ / d'^Xyf^R{g) + J d'^x^/^'R.{f ) 


where and are two dynamical metrics, and R{g) 
and TZ{f) are their Ricci scalars, respectively. The 
parameters Kg = SttG and k^ = 8ttQ are the corre¬ 
sponding gravitational constants, while k is defined by 
K^ = K^ + K^. We assume that the matter action is 
divided into two parts: 


= 5[“l(5,^g) + 4'"'(/,^/), (2.2) 


i.e., matter fields ipg and iljf are coupled only to the g- 
metric and to the /-metric, respectively. This restriction 
guarantees the weak equivalence principle as well as the 
ghost-free condition. The g-matter ipg and the /-matter 
are coupled gravitationally only through the interac¬ 
tion between two metrics g and /. We call ipg and ipf 
twin matter fluids 52|. 

The ghost-free interaction term between the two met¬ 
rics, often called the dRGT potential, is given by 


(2.3) 

'%(7) = 

= (2.4) 

'^ 3 ( 7 ) =-^ We“'''^"7^7V% , 


where are coupling constants, while 7 ^^, is defined by 


7 ^ 7 ^ = • (2.5) 

Taking the variation of the action with respect to g^i, 
and we find two sets of the Einstein equations: 

G'', = + (2.6) 

(2.7) 


where G^i, and are the Einstein tensors for g^^ and 
fgi,, respectively. The matter energy-momentum tensors 
and are given by the variation of matter 

actions, and the 7 -“energy-momentum” tensors 
and are obtained by the variation of the dRGT 

potential action with respect to g^i, and respectively, 
taking the form 25|, 


Th]p^ = ^ 

( 2 . 8 ) 

where = XlLi (- 1 )"“^^ ELo bu+k^^k- 

The energy-momenta of matter fields are assumed to 
be conserved individually as 


( 9 ) . , (/) . , 

= 0 , = 0 , ( 2 . 9 ) 


( 9 ) if) 

where V ^ and V g are covariant derivatives with respect 
to ggi, and fg^. From the contracted Bianchi identities 
for (j2.6l) and (EH), the conservation of the 7 -“energy- 
momenta” is also guaranteed as 

( 9 ) , , (/) . , 

= 0 , = 0 . ( 2 . 10 ) 

The Einstein equations (12.611 and (12.71) , together with the 
continuity equations (12.91) and ( 12 . 1011 . determine the dy¬ 
namics of the bigravity system given in (12.1(1 . In particu¬ 
lar, (12.101) are absent in GR, and in the massive/bigravity 
theory they give non-trivial constraints on cosmological 
solutions. 


B. Homothetic solution 

First we give one simple set of solutions, in which we 
assume that two metrics are proportional; 

fpu=K^gp.: ( 2 . 11 ) 

where A' is a scalar function. This ansatz provides a 
subset of the complete set of solutions, and from the 
energy-momentum conservation (|2.10p . we find that K 
is a constant. As a result, we find two sets of the Ein¬ 
stein equations with cosmological constants Kg and A/: 

Gg,{g) + Kggg, = ( 2 . 12 ) 

Ggu{f)+Kffg, = (2.13) 
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where 

^2 

Kg{K) = (&o + SbiK + 3b2K^ + b^K^) , 

Kj 

^2 

A/(A") = (64 + 363X"^ + 362A:-2 + 61AT-^) . 

(2.14) 

Since two metrics are proportional, we have the con¬ 
straints on the cosmological constants and matter fields 
as 


Kg{K) = K^AfiK), (2.15) 

^ . (2.16) 

The quartic equation (12.151) fixes the value of K. It gives 
at most four real roots, each of which gives a differ¬ 
ent cosmological constant. The basic equations (12.121) 
(or (12.131) 1 are just the Einstein equations in GR with a 
cosmological constant. Hence any solutions in GR with 
a cosmological constant are always the solutions in the 
present bigravity theory. We shall call these solutions ho- 
mothetic solutions because of the proportionality of the 
two metrics. 


III. LINEAR STABILITY ANALYSIS OF A 
HOMOTHETIC SPACETIME 

A. The perturbations around a homothetic solution 

The bigravity theory contains both a massless and 
massive spin-2 fields. This becomes clear when we look 
at the linear perturbations around a homothetic solution. 

The unperturbed solution is assumed to be homoth¬ 
etic, i.e.. 


4 . = , (3.1) 

where bar indicates background quantities. This provides 
solutions to the unperturbed part of the two Einstein 
equations (I2.12|) and (12.131) . A constant K is determined 
by the quartic equation (12.151) . and the matter energy- 
momenta satisfy the unperturbed part of (12.161) . 

We then consider the following perturbations: 

g/iv = g^iv + (3-2) 

= 7m. + (sm. + ^1/i) (3-3) 

where | |, |hj/j | <C suffixes of hjfi as well 

as hjiil are raised and lowered by the background metric 

gfiv 

We obtain the quadratic action for the perturbations 


of the metrics as, disregarding the tadpole-like term, 


82 = 


— £eH Ag] -I-^£eH Ag] 


+ ^\rnls\ 

K _ 

= ^ [ d'^xy/A^CEii[h^~^^;Ag] 

+ ^ J £eh [^^“7 Ag] -b Afp ; m-eff] : 

(3.4) 

where the massless and massive graviton modes are de¬ 
fined by 

- ^{/i . (3-5) 

2 2 

hl+] _ h\g] , fq Q) 


and = Ar“^KgK^/K?_ = k'^ + AT'^Kj are the effec¬ 
tive gravitational constants for the massless and massive 
graviton modes, respectively. The effective graviton mass 
in a homothetic background spacetime is defined by 


:= ^ 



{biK + 2b2K^ + bsK^). (3.7) 


The quadratic Einstein-Hilbert Lagrangian and the FP 
mass term for a metric perturbation are defined by 

CEH[h;Ag] = - T , 

(3.8) 

CEp[h-,mlg] = - ^ , (3-9) 

where 

£liv,aph d = —-dhgi, — m) 

- 2 (gRu)a - -jhRg^ 

-^{gfiuh-2hg^)R. (3.10) 

In the present paper, we assume the graviton mass is 
independent of the cosmological constant, although both 
of them are fixed by giving the coupling constants {bk} 
{k = 0.. .4) and the gravitational constants Kg and k^. 
We treat the graviton mass and the cosmological constant 
as free parameters. It is now explicit from (j3.4p that the 
bigravity theory contains one massless and one massive 
graviton. 
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B. Massive graviton mode in a curved homothetic 
background spacetime 


At the linear order, the massive graviton mode in the 
bigravity theory is given by the FP theory on curved 
background. We note that our background spacetime is 
dynamical unlike the original FP theory. 

The massive mode h\i} does not have a gauge symme¬ 
try, since the mass term breaks the diffeomorphism. The 
gauge symmetry can be explicitly restored by introducing 
the Stiickelberg fields and tt as 

+ 2V(^A) + SV^V^TT. (3.11) 

The perturbation is invariant under the following 
gauge transformations: 

—>■ , Afi —^ Afi — , (3.12) 

and 

Af,^Af, + y^X, (3.13) 

We can interpret Afj., and tt as tensor, vector, and 
scalar graviton modes, respectively. 

The Einstein-Hilbert action preserves the diffeomor- 
phism invariance, and therefore neither .4^ nor tt appears 
in On the other hand, the FP mass term is 

rewritten by using the Stiickelberg fields as 

£fp = - - n ^) - 

o o 

+ ^ - HV^A) 

2 

TTX^ _ _ _ 

(3.14) 

where = 2Vj^A]- Note that the interaction terms 
between the tensor and the scalar graviton modes pro¬ 
duce the van Dam-Veltman-Zakharov (vDVZ) disconti¬ 
nuity around the flat background Q. On the other hand, 
the vDVZ discontinuity does not occur around a curved 
background in the massless limit [s^. 

In order to see the absence of the vDVZ discontinu¬ 
ity, we take a canonical normalization for the tensor and 
vector graviton modes as 




A^ 


Weff 


-An 


and for the scalar graviton mode as 

K- 


^eff \/Rn 


where a positive constant denotes a typical scale of 
the background Ricci tensor, i.e. Rq ^ Then 

the scalar graviton part of the action becomes 


S'2 D 


J 


_ _ R 

,7rV,,7r -t —^ 

A 


TOeff 


R. 


--HOTr) 


(3.15) 


Therefore, when the effective graviton mass is negligi¬ 
ble compared with the background Ricci curvature (i.e. 

<C Rq), the interaction between tensor and scalar 
modes vanishes, and then there is no vDVZ discontinu¬ 
ity. 

However the kinetic term of the scalar graviton is modi¬ 
fied from the standard one in such a massless limit. Hence 
the ghost instability or the gradient instability may ap¬ 
pear, depending on the background Ricci curvature. For 
instance, we consider the background spacetime is given 
by the flat FLRW universe: 


ds^g = a^{—drf + dx^) . 


(3.16) 


The kinetic term is expressed by 
3R2 




2 a 2 


(1 -I- 3w) 


- Y 


w — 1 


3w 




(3.17) 


where w is the effective equation-of-state parameter of 
the background universe defined later by (j4.7ll . and H = 
d/a?' is the Hubble parameter where a dot is the deriva¬ 
tive with respect to the conformal time rj. Therefore the 
Higuchi ghost type instability appears for w < —1/3, 
while the gradient instability is found for —l/3<w < 1 
(^ . This fact indicates that an instability is unavoid¬ 
able if the background spacetime consists of an ordinary 
matter components {w < 1). 

It was shown that the gradient instability appears in 
a more general non-homothetic but homogeneous and 
isotropic background in bigravity 

The instabilities indicate that the linear approximation 
is no longer valid in such a background universe. Al¬ 
though the interaction between tensor and scalar modes 

is suppressed by the small coefficient nies / the in¬ 

teraction between the tensor mode and the scalar mode 
cannot be ignored because the unstable scalar mode 
grows exponentially in time. As a result, unless this in¬ 
stability is resolved, the fifth force will be recovered, and 
then the vDVZ discontinuity will reappear. 

Especially when we discuss the early stage of the Uni¬ 
verse, in which we hope to recover a standard big bang 
universe, this instability becomes a serious problem. If 
the instability and vDVZ discontinuity are not resolved. 
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the bigravity theory cannot realize the realistic homoge¬ 
neous and isotropic universe without an elaborate fine- 
tuned initial condition. 

The above argument is only based on the linear theory. 
The existence of the instability indicates that nonlinear 
interactions cannot be ignored This instability is sup¬ 
posed to arise from the scalar graviton mode. Hence, the 
nonlinear interactions of the scalar graviton mode must 
be taken into account. If the above linear instability is 
stabilized and then the small coefficient proportional to 
the graviton mass term is kept to be small enough, the 
standard big bang universe can be recovered as a stable 
solution in the early stage of the Universe. It is similar to 
the Vainshtein screening mechanism in which the nonlin¬ 
ear effects of the scalar graviton are essential. We shall 
call our case the cosmological Vainshtein mechanism. 

For this reason, we consider the perturbation around 
the flat FLRW background retaining nonlinearities of 
the Stiickelberg fields. We then discuss whether the 
Stiickelberg fields can be stabilized by nonlinear inter¬ 
actions of the scalar graviton, and whether the fifth force 
can be screened in the early stage of the Universe. 


IV. SCALAR GRAVITON WITH NONLINEAR 
EFFECTS 


with 


Ag=K^Af, (4.5) 

<P9 = , (4.6) 

where H = aja} and a dot is the derivative with respect 
to the conformal time rj. We define the effective equation- 
of-state parameter w by 


'^gPg ^g 
^gPg + 


- 1 - 


2H 

SaH^ ■ 


(4.7) 


For general non-linear perturbations, it is difficult, if 
not impossible, to do an analysis even for spherically sym¬ 
metric system without numerical simulations. Hence we 
discuss the following approximated model. First we im¬ 
pose spherical symmetry at full order and assume that 
the g- and /-spacetimes are approximated by the FLRW 
metric such that 


dsl = a^rjg) + e^^^dr^ + r^dfl^] , (4.8) 

ds} = K^a^ipf) [-e^^fdgj + idr) + r)dn^] , (4.9) 

where we introduce two coordinate systems (pg^rg) and 
{pf,rf) to describe the approximated FLRW spacetimes, 
which are given by different coordinate transformations 
from the original one coordinate system ( 77 , r) as 


A. Strategy 

To ease the difficulty in analyzing non-linear effects 
for a generic spacetime, we restrict our analysis to a 
spherically symmetric configuration of cosmological solu¬ 
tions. Even in a spherically symmetric system, however, 
it is still difficult to discuss full nonlinear effects without 
resorting to numerical simulations. Hence, we consider 
some simplified case. In this subsection, we summarize 
the strategy for our analysis. 

We consider non-linear perturbations on homothetic 
flat FLRW backgrounds: 

ds) = a'^{p){-dp‘^+dr‘^+r'^dft^), (4.1) 

ds) = K^a^{p){—dp^ + dr^ + r^did^). (4.2) 

where = dO'^ + sin^ Odip'^. This homothetic solution 
satisfies 


3iL — KgPg Ag, 

(4.3) 

^+3i7(Pg+Pg) = 0 , 
a 

(4.4) 


^ Some nonlinear interactions by introducing the Stiickelberg 
fields on the curved background are discussed in [ 5 ^ . 


Vg=Vgiv,r), rg=rg{p,r), 

Vf=PfiP,r), rf=rf{p,r). 

The approximated FLRW spacetimes mean that we as¬ 
sume |4>g|, likgl <C 1 and |$/|, |4>/| <C 1 because the mass 
interaction term, which is proportional to rri^^g/R q and 
gives the deviation from GR, is assumed be small. How¬ 
ever they do not mean (rig,rg) ~ (? 7 ,r) and {pf,rf) ^ 
(rj, r), in which case the deviation from homothetic space- 
times is small and then can be described by the linear 
perturbations. 

Although the bigravity theory allows one coordinate 
transformation, two independent coordinate transfor¬ 
mations can be possible apparently by introducing the 
Stiickelberg field such that 

+ rf=rgPA'', (4.10) 

where (.4'',.4'') is the Stiickelberg field in the spherically 
symmetric case. Using a gauge freedom, we can fix one 
coordinate system. 

We also assume that for the unperturbed FLRW space- 
times, K and a are given by Eqs. (HU), (g^ and (ITO 
although they are not homothetic. This is allowed as we 
obtain the consistent perturbation equations with this 
ansatz. 

For the following discussions, we consider only a sub¬ 
horizon scale (aL <C H~^) and a length smaller than 
the Compton wave length of the massive graviton (aL <C 
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m~g). This is because our interest is the sub-horizon 
physics during the epoch H > meff and such a solution 
provides us a stable cosmological Vainshtein mechanism. 
We define a dimensionless parameter as 

e:=-^ = aLH, (4.11) 

which satisfies e <C 1 for a sub-horizon scale. 

For a spherically symmetric spacetime, we can divide 
the behaviour of all variables into two: one is a slowly 
changing longitudinal mode mainly due to the back¬ 
ground expansion and matter distributions, and the other 
is a fast changing wave-like oscillation mode of a scalar 
graviton. If the wave amplitude of the oscillation mode 
is small, when we take an average over the typical scale 
of the system oL, which is smaller than the horizon scale 
we find only longitudinal-mode variables, which we 
call adiabatic modes. We then decompose all variables 
X into adiabatic and oscillation modes as 

X = (4.12) 

with 

« (X) , (4.13) 

where ( ) denotes an average over the typical scale of the 
system. 

The dynamical time scale of the adiabatic mode is as¬ 
sumed to be the Hubble time scale, and then its evolu¬ 
tion is caused by the expansion of the Universe and of 
the density perturbations. On the other hand, the os¬ 
cillation mode comes from the degree of freedom of the 
scalar graviton. The time scale of the oscillation mode 
may be the same order of the inhomogeneity scale. Then 
each change rate is evaluated as 

IdrjX^^^l - , (4.14) 

(4.15) 

Since we consider a sub-horizon scale {aL < H~^), the 
dynamical time scale of an oscillation mode is much 
shorter than the Hubble expansion time, i.e., 

» \aHxn ■ (4.16) 

Since the adiabatic mode may be obtained by taking 
the spatial average over the typical scale of the system, 
if the oscillation mode is small, we can assume that the 
dynamics of the adiabatic mode is decoupled from the dy¬ 
namics of the oscillation mode. This assumption is valid 
if the small oscillation mode has no instability. Hence we 
first consider the evolution of the adiabatic modes with¬ 
out the oscillation modes. Then we study the dynamics 
of the oscillation modes around this adiabatic solution. 


B. Adiabatic mode solution 

In this subsection we discuss the time evolution of 
the adiabatic modes for the case without matter per¬ 
turbations in order to see the behaviour of non-linear 
Stiickelberg field. The full analysis including matter per¬ 
turbations will be discussed in Sec. El and the explicit 
expressions will be summarized in Sec IV Al 

For the adiabatic modes, we fix the gauge freedom 
(14.1011 by setting 

119 = 11 , i'g=r ^ ( 4 . 17 ) 

and introduce the dimensionless variables u and /i to 
parametrize 77 / and ry as 

rjf = {I + i,)r], rf = il + fj.)r. (4.18) 

We assume that the time coordinate 77 / and the radial 
coordinate rj point the same directions of rjg and Vg, 
respectively, i.e., 

77>-l, (4.19) 

We have assumed the weak inhomogeneous gravitational 
fields around the FLRW spacetimes, i.e., 

I^9U4'5|,|$/|,|4'/K 1, (4.20) 

and 

|r$;|,|r4';|,|r$}|,|r4''^| <1, (4.21) 

which means that the perturbations from homogeneous 
and isotropic spacetimes are small. Note that this does 
not imply that the perturbations from the homothetic 
FLRW spacetime are small because of the existence of 
the non-linear Stiickelberg field, i.e., either v ot: ji are 
not necessarily small. The nonlinearities in the variables 
v and /i must be retained. However we assume those 
variables are not so large such that the perturbations of 
gravitational fields are still small, i.e., 

|M$g|«l, |V4>,|«1,... . (4.22) 

The spatial derivative of the adiabatic mode may be eval¬ 
uated by 

\drX^’^\ L-^\X^’^\, L<r, (4.23) 

which leads with (14.1411 to 

IdgX^’^l ~ . (4.24) 

Since the dynamical time scale of the adiabatic mode 
variables is given by our spherically symmetric so¬ 

lution around the FLRW spacetimes must be restored to 
the static solution in the limit oi H —0, which is shown 
in Appendix (see also (H). 









In a spherically symmetric static solution, the metric 
f^i, has no non-diagonal component in the coordinates 
{r],r). Hence the non-diagonal component of is at 
most of the order of e in the present adiabatic solution 
around the FLRW spacetimes. The non-diagonal com¬ 
ponent is given by 

frjr = —-I- Tjvyriv' — {r + r^yrfi] . 

(4.25) 

where a prime denotes the derivative with respect to r. 
Because z> ~ aHv, ~ aH^ and |$/|, -C 1, we find 
the leading contribution as, assuming K ^ 0(1), 

fr^r ^ 0{r]iy',efi, ef/) (4.26) 

which must be 0(e). Since rjv' ^ vjt ~ 0(e), the 
Stiickelberg variables v and \i are evaluated as 

IH<0(e^), |mI<0(1)- (4.27) 

We expand all basic equations up to the second order of 
e. The inhomogeneous gravitational fields 
are determined by the Einstein equations, whose explicit 
solutions are given by Eqs. (I5.2|) - (I5.5D in Sec. IV A1 The 
Stiickelberg variable v is solved by the r\ component of 
the interaction conservation law = 0 as 


_9 

dr 


2 -H (r/r) 


- {^v' -I- arHfj.(2 + {rfj.)') + rfj^ 

X (1 + 2/32/i + /33/i") = 0, (4.28) 


where the parameters /32 and (3^ are defined by 


/32 := 


/33 := 


b2K'^ + 

biK + 2b2K^ + baK^ ’ 
b^K^ 

biK + 2b2K2 + bsK^ ' 


(4.29) 

(4.30) 


Eq. ()4.28[l is integrable. An integral constant must be 
zero because of the regularity condition at r = 0. Hence 
we obtain two cases: The first parenthesis in the square 
brackets of (14.281) vanishes or the second one does so. If 


the second parenthesis vanishes, /i is a constant. How¬ 
ever, barring special tuning of model parameters, such a 
solution cannot reproduce the static result in the limit 
of H —0. Hence we conclude that the first parenthesis 
vanishes, i.e., 

r/ji' = —Harfi{2 + (r/i)') — r/i , (4-31) 

which determines v by giving This expression shows 
that the condition (14.271) is consistent. 

Substituting (I4.3ip . together with later obtained 
Eqs. (IO)-(l531). into the r component of = 0, 

we obtain an algebraic equation for another Stiickelberg 
variable ^ as 


Cm 2 (/^) + C //2 (/r) — 0 , (4.32) 

where both Cj „2 and 0^2 are quintic functions of fj, (The 
explicit forms are given in Sec. IV Al) . These terms have 
typical magnitudes given by 

Cm2 ~ miff X , Ch2 - X 0{fi) . 

The equation (14.321) reproduces the static result (IA13|) in 
the limit of —>■ 0. 

Since // is determined by the algebraic equation (14.321) , 
has no dynamical degree of freedom. It is not surpris¬ 
ing because we have ignored the oscillation mode which 
corresponds to the dynamical degree of freedom of the 
scalar graviton. As a result, the Stiickelberg fields do 
not have any dynamical freedom in the adiabatic mode 
solutions. 

From now on, we focus on the period of the Universe 
with H ^ meff, which corresponds to the early stage of 
the Universe. The algebraic equation (14.321) reduces 

Ch2~Q. (4.33) 

This equation has at most four roots, which are given by 
/i = — 1 and 


fiQ = 0, and /r± , (4.34) 


where 


fJ-± = 


1 -I- (I — ■iw)P2 ± V^I — 4/32 + (1 — 3w)2/3| -I- 3(1 — r(;)(I -I- 3ty)/33 

-2/32 + (1 + 3n;)/33 


(4.35) 


r 


Since the root /i = — 1 gives Vf = 0 for any r, we do not 
adopt this solution, and consider only the other three 
roots /xq = 0, /x±. Since those roots are constants, which 


depend on the coupling constants and equation-of-state 
parameter, we can classify the solutions of Eq. (14.321) by 
those roots, which we call the /xo-branches. 
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When we do not include matter perturbations, neglect¬ 
ing the contributions from the interaction terms (which 
are much smaller than for H ^ m-eff), the metric per¬ 
turbations are solved as 


subsection. When the oscillation-mode perturbations are 
defined by 

^osc = X - , (4.41) 


(4.36) 

Then two metrics are given by 

ds^ = a^{r]) [-d?7^ -I- -f , 

dsj = K'^a^{rjf) [—drj'j + dr^ + , (4.37) 

where the coordinates (??/,?'/), which correspond to the 
nonlinear modulations due to the adiabatic Stiickelberg 
fields, deviating from the physical coordinates as 

Vf =r]-^Har^{2fio +dl): rf = {l + fio)r, (4.38) 

where we have integrated Eq. (14.311) for v setting ^ = /ig. 

The solution /rg = 0 corresponds to the homothetic 
FLRW spacetimes, while we can also find the other ap¬ 
proximately homogeneous and isotropic solutions with 
Mo = M± in til® massless limit, in which the coordinate 
transformation from {rj^r) to (?7 /j^/) is nonlinear. The 
solutions with /ig = fj,± are valid up to the second order of 
e with e <C 1 when the interaction terms can be ignored. 
Thus the solutions with the nonlinear Stiickelberg vari¬ 
able M are not exactly homogeneous and isotropic space- 
times, but approximate homogeneity and isotropy still 
hold in the sub-horizon scales. 


C. Stability conditions of scalar graviton 


Next we consider the oscillation modes of perturba¬ 
tions. In the previous subsection, we discussed a spheri¬ 
cally symmetric solution based on the adiabatic potential 
approximation. It does not contain the dynamical degree 
of freedom of the scalar graviton. In this subsection, we 
analyze the stability of the solution against the fluctua¬ 
tions of the scalar graviton with the condition H ^ nieff. 
The stability for the case with H <g; m-eff is discussed in 
Appendix I A 31 

We consider the following perturbations: 


ds^g = a^idg) -I- 

ds} = K^a^ivf) 


+ rjdn^ 


(4.39) 

(4.40) 


We divide the perturbations into the adiabatic and oscil¬ 
lation modes. When we take an average over the typical 
scale of the system, the oscillation-mode perturbations 
do not contribute, and the equations for the adiabatic 
equations are obtained. We solved them in the previous 


the equations that govern their evolution are found by 
subtraction of the adiabatic modes from the full pertur¬ 
bation equations. 

Using a gauge freedom of the oscillation-mode pertur¬ 
bations, as in (I4.10|) . we set two coordinates as 

rig = r] + Sri{ri,r), Vg = r + Sr{r],r), (4.42) 

Vf=V-^Har^{2iio + lil), rf = {1 + iio)r, (4.43) 

where we have used the previous solutions for the adia¬ 
batic mode. 

While ($g, tkg, $/, tk/) are the adiabatic modes, 
{(j)g ,'ipg,(j)f,ijjf,Sr],Sr) are the oscillation modes of pertur¬ 
bations. We assume that all oscillation-mode variables 
have small amplitudes, i.e., ^ 1; a-nd the rate of 

their change in time is roughly ~ |i9i.x°’’'^|- 

We find that the perturbed metric variables are not 
dynamical and they vanish in the limit of rries/H — > 0. 
This is easy to see from the equation of motion as fol¬ 
lows: The Einstein curvature tensors contain the terms 
proportional to H^, while the energy-momentum tensors 
of the interaction term are proportional to rn^ff. For in¬ 
stance, the (Mg, ? 7 g)“Component of the Einstein equations 
in the coordinates (m, r) gives 


6H^g 


2 

t\ _ 


2 d{riljg) 2H dtjjg 

q2j,2 q Qjj 

(1 -f 2/32/io + feMo) djr'^Sr) 

r2 Qj, 


(4.44) 


Taking into account other components of the Einstein 
equations, we see that the energy-momentum tensors of 
the interaction term are negligible compared to the Ein¬ 
stein tensors when 




K_ 




K_ 




(4.45) 


where da = {dg,dr). When the conditions (14.451) hold, 
i.e., for the early stage of the Universe with rrias <C H, 
it is justified that the Einstein equations in bigravity is 
restored to the GR form. In such a stage, the Einstein 
equations for give the solution: 

^g^O, V’g ~ 0 . (4.46) 

This result is convincing because there is no dynami¬ 
cal degree of freedom in a spherically symmetric system 
without matter perturbations in GR. By the same argu¬ 
ment as above, we also find 


V'/^O, (4.47) 
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from the Einstein equations for /^i,. 

In the limit of mes/H —>■ 0, we find that two metrics 
without matter perturbations are given by 

dsl = a^{r]g) [-d'q'^g + drl + rldn^] , (4.48) 

dsj = K^a^{rif) [—dr]j + dr^ + r^dfl^] . (4.49) 


We then expand the action in terms of {Srj, dr) up to the 
second order of e. The variations with respect to Sr] and 
Sr give the constraint equations, which are solved such 
that the Stiickelberg variables Srj and Sr are given by 


Sr] = 


Sr = 


1 + Mo 

^r-TT arH]lo dr,TT 2', 

a2(l+/io) 1 + Mo a2 


(4.50) 

(4.51) 


in terms of a Stiickelberg scalar tt. For the analysis of 
stability, it is sufficient to determine Srj and Sr up to the 
first order of e. 

Substituting (j4.50l) and (14.511) into the action, we ob¬ 
tain the quadratic action of tt as 


& = 


'■eff 


dfl 


dr]dr{arH)^]Cs (9^7r)^ — Cg {drir)^ , 

(4.52) 


where fl is the solid angle. The signs of these coefficients 
K-s and c| determine the stability of the Stiickelberg 
scalar tt, which corresponds to the scalar graviton. 

For the mo = 0 - branch, the coefficients are given by 

»Cs|„.o = j(l + 3»i), (4.53) 

which is consistent with the result (13.171) . On the other 
hand, for the mo = M±" branches, after simplifying the 
expressions by using (14.351) . we find 

^sUo=Ai± = ~ 1)(2 + M±)(l + 2/32M± + /33M±) 1 

(4.55) 

2 I ^ 2 (3(1 - w) -b (1 - (3w - 1)/32)m±) 

CsIa4o=a4± - 3(3^ _ + ^±)(1 + 2/32M± + ' 

(4.56) 


The no-ghost condition is given by /Cg > 0 while the no¬ 
gradient instability condition is given by Cg > 0. Hence, 
the stability condition of the scalar graviton is to have 

/Cg > 0, c| > 0 . (4.57) 


In the Mo = 0-branch, either the ghost instability or 
the gradient instability appears for w < 1. So the in¬ 
stability is inevitable when the Universe consists of the 


standard matter, as shown in Sec. IIIIBI Since the insta¬ 
bility appears only in the relation between the two coor¬ 
dinate systems, (MgUg) and ( 77 /, r/), two spacetimes still 
keep homogeneous approximately as long as the condi¬ 
tion (14.451) holds. However, since tt grows in time due to 
the instability in this branch, the condition (14.451) even¬ 
tually breaks down. 

On the other hand, the mo = M±" branches can avoid 
the ghost instability as well as the gradient instability de¬ 
pending on the background dynamics and the coupling 
constants. In Fig. (TJ we show the parameter regions 
where the solution is stable for the cases of w = — 1, 
w = 0, w = I/3± 10“®, w = 2/3 and w = 1 —10“®. Note 
that in the radiation dominant universe with w = 1/3, 
the action is given by 

‘^ 2 L=i /3 = J J dr]dr{arH)^ {drirf , 

(4.58) 


which does not describe the dynamics of tt Thus the 
expansions we have adopted in our calculation is inval¬ 
idated in this limit, and in order to correctly study the 
dynamics of the scalar graviton, we must calculate the 
higher-order terms of tt. 

For w ~ 1/3, the existence of the stable solution is 
guaranteed for the parameter region such that 

I3l> /33> 132, for w = ^- |(5w|, (4.59) 

Ids > (32, for w = i -I- |(5w|, (4.60) 

with |(5r(;| <C 1. In such a parameter region, at least one 

of ]i± satisfies the stability condition (14.571) as well as our 
ansatz ]j,± > — 1. 

For Tc = 1, ]i± are given by 


M± 


1 - 2/32 ±|1-2/^21 
2(2/33 - (32) 


(4.61) 


Hence, one of m± becomes zero, which gives the homoth- 
etic solution. 

Since reducing the stability condition (j4.571) to the al¬ 
lowed parameter region for arbitrary values of ic is a non¬ 
trivial task due to the complicated dependence (14.351) of 
M± on the model parameters, we analyze the stable re¬ 
gion numerically. We conclude that the parameter region 
of 


(3l >(33 > 132, (4.62) 


If we take into account the trace anomaly of quantum cor¬ 
rections, we find small deviation from w = 1/3. For example, 
_ 1 5 .<,* (JV„ + fjV^)(JijVe-§JV^) 

SU(Afc) gauge theory with coupling g and Nf flavors ISfiH . 
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FIG. 1: The stable regions of the coupling constants for flat FLRW backgrounds for the cases of w = —1, w = 0, w = 1/3±10“®, 
w = 2/3 and w = 1 — 10~®. The fj,o = fJ.+ - branch is stable in the regions denoted by fj.+ (red regions), the /ro = “ branch is 

stable in the regions denoted by fi- (blue regions), and both fio = fj,± are stable in the regions denoted by fj,± (purple regions). 
The black solid curves correspond to ps = P 2 and the black dotted lines correspond to P 2 = Ps- 


guarantees the existence of a stable branch for any values 
of w except for w = 1/3. Even outside the region (14.621) . 
we obtain stable branches for some values of w, but the 
instability always appears for near radiation dominant 
stage such that w = 1/3 — |(5w|. 

In the stable region (14.621) . only one stable branch ex¬ 
ists for any w. When the Universe consists of the usual 
matter field with w < 1 , either /tq = Ms- or /i_ gives a 
stable solution in the parameter region (14.621) . While if 
the Universe is composed effectively of a “strange” mat¬ 
ter field with w > 1 , we find ^0 = ft±-branches are 
stable only in the region of P 3 > /3| (see Fig. El, and the 
stable region disappears in the limit of w = 00 as shown 
in Fig. El As a result, only the /tq = 0 -branch becomes 
stable. 

Although the condition (j4.62() depends on the propor¬ 
tional factor K, we can derive the stability condition for 
the original coupling constants {bi}. Since we assume 
TOgjj > 0 (i.e. Kbi + 2 K %2 + K^bs > 0 ), the condition 




ft ft 

(g) ■!« = H- 10“^ (h) w = 2 

FIG. 2: The same figures as Fig. [T] for the cases of w = 
1 -I- 10“® and w = 2. 


reduces to 

-6163 >0, 62 <0. (4.63) 

Therefore, if we choose the coupling constants in the 
above parameter regions ()4.63l) . both ghost instability 
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and gradient instability are avoided because of the non¬ 
linear interactions, and then the early stage of the Uni¬ 
verse described approximately by GR FLRW spacetime 
becomes stable. 

V. PERTURBATIONS WITH MATTER 
EFFECTS 

In this section we consider the adiabatic modes re¬ 
trieving matter perturbations and discuss the evolution 
history with the nonlinear effects of the Stiickelberg 
fields. For the stable branch discussed in the previous 
section, the evolutions of the matter perturbations are 
approximated by the adiabatic modes. The oscillation 
modes of matter fluctuations decay in time [s^. Hence 
we can decouple the adiabatic and oscillation modes in 


this case as well. 


A. Adiabatic mode with matter perturbations 

Now we discuss the time evolution of the adiabatic 
modes with matter perturbations. The matter energy 
densities are perturbed as 

Ps =Ps(l + <^s): P/= P/(l + <^/)- (5-1) 

We ignore pressure perturbations and spatial velocities 
compared to the density perturbations, just for simplic¬ 
ity. 

All equations are expanded up to the second order in 
e. The Einstein equations for and give 


and 


2^'g(77, r) = a^{r])r^ rnl(g + 




I I 3^ I ^ 2 - r 


2^f(Vf,rf) = a^(vf)rf 


d^f 


2 ’’/y:^(P/U/) = a (gf)rf 


(1 -I- 


(1+Mr 


^ ^ P + (1 + /32)p^ + ^ ^ fJ-A + ^K^^jpfSf 


2 I 2 -|- 2^2 /^3 3'\ 1 J^2 2 - J 

H + 2fl + ---y -t -K KfPfSf 


(5.2) 

(5.3) 

(5.4) 

(5.5) 


respectively, where 


hiv,r) 


Jo 47rf^(5gdf 


^/(P/U/) 


fo^ 47rr^Sfdf 
Airr'^dr 


(5.6) 


are spatial averages of the density perturbations in the spheres with the radii r and rj, respectively. We define the 
mass parameters by 


ml-.= '^^^{biK + 2b2K^ + b 3 K^), m} := + 2b2K^ + bsK^), (5.7) 

hh t\ tv 

with which the effective graviton mass is expressed by and as = nv^ + m^. 

Although the /-variables are given as the functions of (??/,?'/), it is easy to find them as the functions of {g^r) by 
use of the Stiickelberg fields v and g. The variable v is determined by (I4.31|) even when the matter perturbations are 
included. 

Substituting Eqs. (I5.2p - ()5.5I1 and (j4.31|) into the r component of = q, we obtain an algebraic equation for 

fi: 


Cm'^ (/i) -K Ch2 (m) + C matter (p) = 0. 


( 5 . 8 ) 
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where we define each function as 

Cm2(/^) := + nf [9 + 18/?2/r + (6/3| + 

+ TTij 9 + 18(1 + /32)/r + (10 + 34/?2 + d/Jg + 6/3|)/i^ 

+ (2 + ld/32 + 8/3g + \20^)ij^ + (2/32 + 2/3g + 2/3^ — /3g + d/32/3g)/i^ | 
•= —3i?^//(l + {3(1 — w) + 2 [1 + (1 — 3 ?c)/ 32] // + [2/32 ~ (1 + 3rc)/3g] /i^} , 

Cmatter(lt) := (1 + {«^g/Og4(l “ -^3/^^) “ K^KjpfSf{l + p) [l + 2p + (2^2 “ | • 


r 


(5.9) 

(5.10) 

(5.11) 


Eq. (15.8t) reproduces the static result (IA13I) in the limit 
of i? ^ 0. 

An important difference from the case without matter 
perturbations is that the variable p depends also on the 
matter perturbations. When the matter perturbations 
are not negligible, Eq. (Id.32p should be replaced with 
Eq. (15.8L in which the typical value of the additional 
third term is evaluated by 

^matter ^ ^gPg^g ^fPf^f • 

The metric perturbations are given by the GR results 
with the corrections coming from the interaction term, 
e.g., one of the perturbations of g^i, is given by 


First, we consider the case of w < 1. The stable branch 
is given by one of p±, and thus (15.131) is also solved as 
P ~ P±- The gravitational sector is restored to the one 
in GR, when 


— < 5 


(5.15) 


i.e., if the correction terms from the graviton mass are 
negligible compared to the GR terms in Eqs. (lOl-dS^I). 
As a result, in the early stage of the Universe, the metric 
perturbations are restored to the GR results. 

Next, we consider the case of ic > 1, in which the stable 
branch is given by /tq = 0. In this case, the solution is 
given by, assuming |/i| -C 1, 


Tg = 4 'gr + X ^ X C>(/i). (5.12) 

Kj _ 

When the second term is negligible compared to the first 
one, the metric perturbations are restored to the GR re¬ 
sults. Since the equations of motion of twin matters are 
not modified from usual ones (e.g., see Eq. (12.91) 1. the 
restoration of the metric perturbations guarantees the 
dynamics of the matter is also restored to the GR result. 
Therefore we will discuss only the metric perturbations 
which are determined by p, as in Eqs. (IO)-(l531). 


B. GR phase 


^g 

3(1 — w) ’ 


(5.16) 


where we have used the background equations. One of 
the metric perturbations is described as 


2-^g = a^r^H^Sg 


2 2 2 


^9 h-h 

3(1 — w) 


(5.17) 


Since the second term is negligible compared to the first 
term in the case oi H ^ rrieg, the metric perturbations 
are again restored to the GR results for w > 1. 

Hence both cases show the GR limit in the early stage 
of the Universe {H ^ rries). We shall call this stage the 
GR phase. 


We discuss a stage when the Hubble parameter is larger 
than the effective graviton mass (i.e., ^ 

cause Cjj 2 Cm 2 , Eq. {El becomes 

Ci/2+Cmatter~0. (5.13) 

Since the second term is much smaller than the first term, 
Eq. (15.131) is schematically solved as 

p = Po + O{5gjf), (5.Id) 

where |(5g|, |(5/| 1. As discussed in Sec. Hvl the stable 

branch is found with po = p± for w < 1, or with po = 0 
for w > 1 . 


C. Bigravity phase 

Secondly, we discuss the stage when the Hubble pa¬ 
rameter is smaller than the effective graviton mass (i.e., 
<C TOgjj). In this stage, we find 6^2 Ch^- Hence, 
for the matter of our interest, Eq. (15.81) reduces to 

+ ^matter ~ 0 . (5.18) 

We denote the roots of Cm^{p) = 0 by poo, which are 
found to be zero and some constants of order unity. Sim¬ 
ilarly to the previous subsection, Eq. (j5.18l) is solved as 

M — Moo “1“ 2 ^ ^(^57 ^/) * 


(5.19) 



















14 


Since Eq. (15.181) is a polynomial equation of degree 
seven for /i, there are seven solutions for fioc- We classify 
the solutions of Eq. (15.181) into two types: the linear 
branch and non-linear branches. Note that a branch here 
denotes one with fi = fioo in the limit of H/rUeS —t 0 . 

The linear branch is realized by choosing Hac = 0. Eq. 
(15.181) gives the value of /x as, assuming |/i| <C 1, 

__ 4 p,s, - lApfS, 

9m,y ■ 

Substituting this solution into the expression of the grav¬ 
itational force (lOl) . we find 




a^r 



^gPg 


3mlg 


K‘ 


KfPf 


(5.21) 


Hence the gravitational force in the g-sector is produced 
by the /-matter as well as the g-matter. We have done a 
detailed study including stability about this linear branch 
in our previous paper [30j. 

The nonlinear branches are obtained by choosing pco 
0 ( 1 ). In this case, since |/ioo| ^ |5g|,|5/|, the solutions 
are found to be 


P ~ Poo 1 (5.22) 

giving the same as those without matter perturbations. 
For these nonlinear branches, the metric perturbations 
include the correction terms of the graviton masses as 
given in Eqs. (I5.2I) - (I5.5I) . Hence, for these branches, there 
is a non-negligible inhomogeneity at large scale and the 
gravitational behaviour deviate largely from GR’s one 
beyond the Vainshtein radius as shown in Appendix lA 31 
We shall study the stability around the nonlinear 
branches as well as the linear branch in the bigravity 
phase (Wgjj ^ H^). In the limit 77 —>■ 0 or —>■ Ag/3, 

we find the adiabatic solutions turn to be the static ones 
given in Appendix |Al In the sub-horizon scale, the cos¬ 
mological coordinates ( 77 , r) are related to the static co¬ 
ordinate {T,R) as 

T = t +—a^r'^■^J~^ + ■ ■ ■ , R = ar , (5.23) 


where ^ ~ j cosmic time. On the cosmological 

coordinate (77, r), the Stiickelberg variable p, is defined by 


p := 




(5.24) 


while on the static coordinate (T, R), p is defined by 

Rf - R 


p := 


R 


(5.25) 


Because of 


p = p, (5.26) 

for Eq. (I5.23|l , we find Eq. (15.81) corresponds to Eq. (IA13I) 
even in the limit 77^ —>■ Ag/3. 

Therefore, as a lowest-order approximation, we re¬ 
gard the adiabatic solutions as the static ones and we 
can apply the stability analysis in Appendix IA3I to the 
present case. We note that the asymptotically homoth- 
etic branch of Eq. (IA13I) corresponds to the linear branch 
in the bigravity phase, while non-asymptotically homo- 
thetic branches correspond to nonlinear branches in the 
bigravity phase. It turns out that the linear branch is 
always stable. Furthermore, we find that one of the non¬ 
linear branch is stable for some coupling constants. 

These arguments are essentially unchanged even 
when a cosmological constant with the condition of 
Ag < 3m'^ffl2 is introduced. We give the results with such 
a cosmological constant for the gravitational behaviour 
and the Vainshtein screening in Appendix I A 21 and those 
for the stability in Appendix lA 31 

Since the gravitational behaviours are modified from 
the ones in GR due to the existence of the fifth force, me¬ 
diated by the scalar mode of graviton, for both branches, 
we call this stage the bigravity phase. 


VI. TRANSITION FROM GR TO BIGRAVITY 

As we have shown in the previous section, the Uni¬ 
verse in the bigravity theory has some stable branches in 
the both limits of 77 ;g> rrieff and of 77 ^ mes, which 
correspond to the early and late stages of the Universe, 
respectively. More precisely, in the early stage of the 
Universe (77 ^ rUes), the po = 0-branch is unstable, 
while Po = p±- branches can give the stable GR phase 
depending on the coupling constants. On the other hand, 
in the late stage of the Universe (77 <C nies), the branch 
with \p\ ^ 1 , which provide us the bigravity phase, is 
stable for appropriate coupling constants. In addition, 
some branches with p^o 0 ( 1 ) (the bigravity phase) are 
also stable for some coupling constants. 

The question is whether any two stable branches with 
different limits can connect under our adiabatic approx¬ 
imation or not. We shall discuss this possibility in this 
section. 

Since our Universe is homogeneous at large scale, the 
Stiickelberg variables should reach to the linear bigravity 
phase as the Universe expands. Hence our Universe must 
start from the GR phase and transit to the bigravity 
phase. Then the Universe must pass through the period 
of 77 ~ TTXeff, where the behaviour of p becomes unclear. 

One unknown in this period is the transition time scale, 
if transition occurs. If the transition time scale is given 
by the Hubble time scale and the adiabatic mode is an 
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attractor even in the period of i? ~ meff, we can discuss 
the transition by considering only the adiabatic modes. 
However, if two stable adiabatic solutions in GR and in 
bigravity phases are discontinuous, the adiabatic approx¬ 
imation breaks down, and then the transition time scale 
may be faster than H~^ even if such a transition is pos¬ 
sible. As a result, the assumption (14.241) is no longer 
valid, and a full analysis without approximation will be 
required. 

Here, we only speculate some possible transitions based 
on the adiabatic approximation with the assumption 

(14.241) . When the amplitudes of the density perturbations 
are given, one can obtain the Stiickelberg modulation 
in terms of H and Sg/f by solving the algebraic equation 
(15.81) together with the use of (14.31) and (14.61) . Since H de¬ 
creases in time, we find the evolution of fj, without solving 
the equations of motion, under the assumption that the 
solutions to (15.81) at different moments are continuously 
connected. Here we also assume that the density pertur¬ 
bations are constant in time; although they may evolve 
in time, the qualitative behaviour does not change much. 

We show several examples in Figs. [3] and 21 We have 
explored a wide range of the parameters with the stability 
condition (I4.62p . The stability analysis for the late stage 
of the Universe {rriee/H 1) is given in Appendix lA 31 

As shown in Figs. |3]and|4l two stable branches in the 
limits of nies/H <C 1 and of mes/H ^ 1 can be con¬ 
tinuous or discontinuous depending on the equation-of- 
state parameter w and the density perturbations. For 
example, as shown in Fig. [3] (H), for w = 2/3, we find 
one continuous curve from the early stage to the late 
stage of the Universe, exhibiting the branch that is sta¬ 
ble in both limits of niee/H <C 1 and of nies/H ^ 1, 
if 5g{= 10“^) > Sf{= 10“^). On the other hand, if 
Sg{= 10 “^) < Sf{= 10 “^), two (one stable and another 
unstable) continuous curves are splitted into two discon¬ 
tinuous curves as found in Fig. [3] (H). There are points 
at which two real roots of Eq. (15.8|) degenerate and be¬ 
yond which they become imaginary. At such points, the 
time derivative of fi {dfx/dH~^) diverges, but fi itself is 
finite. At such a singular point, we argue that the en¬ 
ergy density of the scalar graviton diverges. This is be¬ 
cause the variable fi is basically related to the adiabatic 
mode of the scalar graviton ttq as ^ = a^^TTp/r. Similarly 
to the Vainshtein screening mechanism, the derivative of 
rfi ddno) gives a typical energy density of the scalar 
graviton. So when dfx —?► ±oo (or ddn —>■ oo), a singular 
behaviour of the scalar graviton may appear. 

However, near such a singular point, our assumption 

(14.241) is no longer valid. This singular behaviour may 
simply be an artifact of our adiabatic assumption, and if 
we analyze the full nonlinear dynamics without the adi¬ 
abatic approximation, this singularity may not appear. 
Hence, when we find a discontinuity in the solution n{H), 
what happens in the transition period is still an open 
problem. In what follows, we just discuss the adiabatic 
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FIG. 3: Case A: There exists the stable branch with non¬ 
linear /ioo(~ 0(1)). We plot the roots of Eq. (15.81) for Sg = 
10~^,5f = 10~^ (red curves) and for Sg — = 10“^ 

(blue dashed curves). We set /32 = —3, ps = 3, rUg = mj. The 
branch with fio ~ 0.40 for ic = 0 and fio ~ —0.10 for w = 2/3 
are stable in the early stage of the Universe (nies/H -C 1). 
In the late stage (mesIH 3> 1), the branches with /tq ~ 0 
and 0.70 are stable. The Universe may evolve from the stable 
/io - branch to the stable ^oo - branch. For example, for w = 0, 
^ changes from po = 0.4 (GR phase) to fiao = 0.7 (nonlinear 
bigravity phase), while for w = 2/3, it does from fiQ = —0.1 
(GR phase) to /Too = 0 (linear bigravity phase). 


solutions given by the continuous curves. 

We classify our results into two cases by whether the 
stable nonlinear branch in the late stage of the Universe 
exists (Case A) or not (Case B). Some examples of Case 
A are given in Fig. O while those of Case B are in Fig. 21 
We shall discuss them in due order. 

(1) Case A: 

For w < 1/3 (Fig. 21(1)), there exists a stable /rp-branch 
in the early stage of the Universe {H ^ mes), and non¬ 
linear stable /Too -branch in the late stage {H <g; nies)- 
Hence we expect that the Universe will evolve from GR 
phase to non-linear bigravity phase. The GR phase in 
the early stage transits to the inhomogeneous universe in 
the late stage. 

Our present stability analysis, however, is valid only 
for both limits ot H ^ rUeS and oi H uieff. There- 
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(11) w = 2/3 

FIG. 4: Case B: There exists only one stable branch with 
/Too = 0. We plot the roots of Eq. (15.81) for Sg = W~^,Sf = 
10~® (red cnrves) and Sg = 10“®,5/ = 10“^ (bine dashed 
curves). We set Pi = —1,^3 = lf2,m'^ = rrif. The stable 
branches for rrigg/H <C 1 are given by fiQ ~ 1.1 for w — 0 , 
and fj .0 ~ —0.21 for w = 2/3. For m^fs/H ^ 1, only — 0 
is stable. The Universe with w = 2/3 may evolve from fio — 
—0.21 (GR phase) to Hoo = 0 (linear bigravity phase), but for 
ui = 0, there is no stable adiabatic solntion. 


fore, it is not immediate to decide whether our adiabatic 
solution is still stable at ~ nies- If the adiabatic solu¬ 
tion becomes unstable near the period of i? ~ rUeS, the 
adiabatic mode does not evolve independently of the os¬ 
cillation mode. Then the transition time scale becomes 
faster than H~^. As a result, the GR phase in the early 
stage might connect to the linear bigravity phase in the 
late stage. 

On the other hand, in the case of w > 1/3, only for 
the case of 6g > Sf, there exists only one stable branch 
from ^0 to /ioo ■ It may describe that the GR phase in the 
early stage of the Universe transits to the linear bigravity 
phase in the late stage. If < 5/, two stable branches 
are disconnected. The assumption (14.241) may no longer 
be valid at the singular points. 

(2) Case B: 

For the case of w > 1/3, we find similar behaviours to 
Gase A. There exists one stable branch from the GR 
phase in the early stage of the Universe transits to the 


linear bigravity phase in the late stage, if (5g > 5/. If 
Sg < 5f, however, two stable branches are disconnected. 

For w < 1/3, the stable GR branch evolves into a 
singularity. No other stable branch exists. Since the as¬ 
sumption (14.241) may no longer be valid, the transition 
from the GR phase to the linear bigravity phase might 
be possible, if the stable /Too = 0 - branch is a global at¬ 
tractor. 

From the above analysis, we conclude that for w > 1/3, 
there exists a stable adiabatic solution, which describes 
that the Universe evolves from a GR phase in the early 
stage to a linear bigravity phase in the late stage. For 
more realistic equation-of-state parameter w < 1/3, how¬ 
ever, if the transition from the GR phase to the linear 
bigravity phase occurs, the transition time scale is likely 
to be faster than H~^ and then the adiabatic condition 
must no longer be valid In order to conclusively an¬ 
alyze the full evolution history of the Universe for nat¬ 
ural equation-of-state parameter re, the analysis beyond 
the adiabatic approximation is required, and a numeri¬ 
cal treatment should be expected, which is the work in 
progress. 


VII. CONCLUDING REMARKS 

Since the homogeneous and isotropic universe is consis¬ 
tent with cosmological observations, the existence of ex¬ 
ponential instability around the FLRW spacetime would 
ruin the viability of such cosmological solutions, unless 
the initial condition is fine-tuned. It has been known that 
in the bigravity theory, while the FLRW background so¬ 
lutions are stable against small perturbations at the late 
stage of the Universe (H m-eff), the linear perturba¬ 
tions suffer either a ghost or a gradient instability dur¬ 
ing the early epoch {H mes). However, the unstable 
mode arises from the Stiickelberg fields which determine 
the relations between the two coordinates of the physi¬ 
cal and fiducial sectors in the bigravity. The instability 
implies that the nonlinear interactions of the Stiickelberg 
fields may be signiheant. However, the nonlinearities of 
the Stiickelberg fields do not directly conclude that the 
spacetime perturbations around the homogeneous and 
isotropic spacetime are large. For this reason, we study 
small metric perturbations around the flat FLRW back¬ 
ground retaining nonlinearities of the Stiickelberg fields 
and discuss whether or not the viable cosmological solu¬ 
tion exists even in H ^ Weff. 


® One may argue tlie possibility that the transition might occur 
with a Hubble time scale by taking into account the evolutions 
of the density perturbations. However this is unlikely because 
the qualitative behaviours of Figs. [3] and |4] do not depend on the 
amplitudes of Sg and Sf but rather on their ratio Sg/5f. 
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In order to comfort the difficulty in the full nonlinear 
analysis, we consider the following simplifications: 

1. We assume both background spacetimes are almost 
isotropic and homogeneous. This background is re¬ 
garded as a ground state of the massive graviton. 

2. Since the scalar graviton mode is essential to the 
Higuchi ghost and the gradient instability, we re¬ 
strict our analysis to spherically symmetric config¬ 
urations, in which the degrees of freedom of the 
tensor and the vector gravitons do not exist, while 
the scalar graviton propagates. 

3. We decompose the perturbation variables into adi¬ 
abatic modes and oscillation modes, and assume 
that the adiabatic modes evolve independently of 
the oscillation modes. The later assumption is valid 
only if the oscillation modes are stable. 

4. We only consider the sub-horizon scale and as¬ 
sume that the adiabatic modes satisfy the condition 
(I4.24L Under these assumptions, we find that the 
result in the adiabatic modes reproduce the static 
result (IA13I) in the limit H ^ 0. 

Then, we show that the Higuchi type ghost instability 
and the gradient instability can be removed by the non¬ 
linear effects of the scalar graviton. In the early stage 
of the Universe {H ^ m-efr), we find a few sets of cos¬ 
mological solutions (I4.37|) . in which there are two types; 
one is given by the homothetic spacetimes, and other 
branches are given by the approximately homogeneous 
and isotropic spacetimes on the different coordinates. In 
the former homothetic solutions, the instability cannot 
be avoided when the Universe consists of ordinary mat¬ 
ter components. On the other hand, the solutions in the 
other branches can avoid both instabilities if the coupling 
constants satisfy the simple inequalities (14.6311 . For any 
values of the equation of state, there exists at least one 
stable branch, in which the Universe evolves in the same 
way as in GR. The result suggests that the bigravity has 
the healthy massless limit, and cosmology based on the 
bigravity is still viable even in the early stage of the Uni¬ 
verse. 

Hence we conjecture the following cosmic scenario: 
First, the cosmic growth history is restored to the GR 
result in the early stage {H ^ m-eff), in which the 
Stiickelberg helds are non-linear. After a transition 
around H ^ meff, the Universe reaches the bigravity 
phase, in which the modification from GR appears. In 
the bigravity phase, the ^oo = 0-branch is stable. There 
is also another stable branch for some combinations of the 
coupling constants, showing somewhat peculiar gravita¬ 
tional behaviours. 

The transition from the GR phase to the bigravity 
phase is, however, still an open problem. One important 
key point is the transition time scale. As discussed in 


Sec. Ell if the GR phase connects to the linear (/Too = 0) 
bigravity phase with natural equation-of-state parame¬ 
ter {w < 1/3), the transition time scale should be faster 
than H~^. To analyze this transition phase, a numerical 
treatment should be required. 

As we have shown, the GR cosmology is recovered in 
the early stage of the Universe, which we call the cosmo¬ 
logical Vainshtein mechanism. Since this is different from 
the Vainshtein screening mechanism discussed in some 
scalar-tensor type gravity theories such as the Galileon 
model, here we show how they are different. 

In the present bigravity model, when we take a decou¬ 
pling limit, we find the effective action of a scalar field </>, 
which corresponds to the massive scalar graviton, as 


= 


-b 


/ 




’(0) 


Ad 




'2 w A.-NL 

^/xt/pcr 


m 


eff 
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eff 


Ad 

“'NL/ii-'pcr 


where 


(7.1) 


= (7.2) 


where cnl,cnl are some dimensionless constants, and A 
is a typical scale for a strong coupling state. We have 
dropped higher interaction terms. In the conventional 
Vainshtein screening, the first two terms play the essen¬ 
tial roles. When the energy scale of the scalar field is 
larger than the string coupling scale A below the Vain¬ 
shtein radius, the degree of freedom of the scalar held is 
killed by the second non-linear term , and then GR is 
recovered. On the other hand, in the present cosmologi¬ 
cal Vainshtein screening, the third and fourth terms play 
important roles. We Hnd the GR solution in the early 
stage of the Universe as follows: Gomparing the hrst and 
third terms, when the typical scale of the curvature is 
larger than the graviton mass, we Hnd the GR universe 
as a consequence. However, to stabilize this GR universe, 
we have to take into account the fourth non-linear term 
as well. Hence in the cosmological Vainshtein 
screening, the curvature terms are essential, which is dif¬ 
ferent from the conventional Vainshtein mechanism. 

Our Vainshtein mechanism is in some respects similar 
to the ghost condensate [57|. In the ghost condensate, 
even if a scalar Held has a ghost instability at linear level, 
non-zero time derivatives of the background expectation 
value can stabilize the fluctuation of the scalar Held. In¬ 
deed, in our system, assuming that the scalar graviton 
(j) can be decomposed into the adiabatic mode ttq and 
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the oscillation mode tt as (/) = ttq + tt, the non-zero ttq 
stabilizes the fluctuation tt. 

However there are some differences between the ghost 
condensate and the stabilization in the bigravity. First, 
we note that the non-zero ttq can stabilize not only the 
ghost instability but also the gradient instability in the 
bigravity. Furthermore, the scalar graviton is stabi¬ 
lized by the higher order non-linear effects of the spatial 
derivative of ttq as distinct from the time derivative in the 
ghost condensate. This is because is related to ttq as 
fiQ = oi~‘^drT^Qlr, and the higher-order terms of are es¬ 
sential to the stabilization. Although the scalar graviton 
has inhomogeneity, the spacetime is still almost homoge¬ 
neous due to the existence of the screening mechanism. 

In order to check complete viability of the cosmolog¬ 
ical model in the bigravity, we must consider beyond 
spherically symmetric perturbations as well as beyond 
our simple ansatz on the background. The kinetic term 
(I3.17|) indicates that the ghost or the gradient instability 
may also exist even for non-spherical small perturbations 
around the FLRW background. Thus extended studies 
are required. 

The stabilities of the vector graviton and the tensor 
graviton are also non-trivial questions. The gravitational 
field equations are restored to the Einstein equations in 
GR when H ^ Weff- Hence, the tensor mode should be 
reduced to the GR one, and then we expect that no in¬ 
stability may appear, just as in GR. We also expect no 
exponential instability in the vector graviton, because the 
kinetic term for the vector graviton is not modified even 
for a curved background at linear level as shown in the 
equation (|3.14l) . The complete stability analysis, how¬ 
ever, should be done to confirm our speculation, which 
we leave for future work. 
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1. Basic equations 


Spherically-symmetric static spacetimes in the bigrav¬ 
ity theory are assumed to take the two metrics of the 
form 


dSg = — e 

dsj = 


+ e^^odRl -f Rldn^ , 
-e^^fdTf + e^^fdR} + Rjdfl 



(Al) 

(A2) 


In what follows, we fix the gauge as 


Tg=T, Rg=R. (A3) 


Since the spacetimes are static, a non-diagonal compo¬ 
nent of on the coordinates (T, R) should vanish to 
find a non-trivial solution ^ Hence Tf must be 
proportional to T. In addition, we expect two space- 
times are asymptotically homothetic at large distances 
(i? —> oo), where the proportional factor K is given by 
Eq. (14.5|) . So we set 


Tf=T, Rf = R[1 + fi{R)] , (A4) 


with the property fl{R) —> 0 in the limit R ^ oo. 

As in [ 5 ^ . we consider the region inside the Comp¬ 
ton wavelength of the massive graviton and weak grav¬ 
itational fields such that the metrics and their spatial 
derivatives satisfy the following conditions: 


|I>,|,|^,|,|I>/|,|'I'/|«I, (A5) 

|R$;i,|i?'F;|,|i?l>'/|,|i?'F'/|«I, (A6) 


Appendix A: Vainshtein screening in a 
spherically-symmetric static spacetime 

In this Appendix, we consider the Vainshtein screen¬ 
ing mechanism in a spherically-symmetric static space- 
time. While the Vainshtein screening in bigravity for a 
spherically-symmetric static spacetime was discussed in 
[55|, the effects due to the existence of the /-matter and 
cosmological constant have not been taken into account. 
We would like to include these effects in our present anal¬ 
ysis. 


where a prime denotes the derivative with respect to R, 
which is used only in this Appendix. Although the vari¬ 
able fl is not necessarily small, we assume that fi should 
not be so large, satisfying 

|A$g|«I, |i?A'^g|«I, •••. (A7) 


From the Einstein equations, we obtain 
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(A 8 ) 

(A9) 

(AlO) 

(All) 


where we introduce the mean densities in the spheres with the radii R and Rf hy 


rRf 


Pg{R) = 


4:TTR^Pg{R)dR 


rR 


pfiRf) = 


10 


AttR^ PfdR 


rRf 


(A12) 


/ AiiR^dR / AnR^dR 

'0 Jo 


with pg = —and Pf = —respectively. We ignore the pressures of twin matters, for simplicity. 

r 1 

Substituting them into V = 0, in the weak field limit, we find an algebraic equation for p as 

^m?^p') “t“l'A(A) “t“ tlniatter(/^) — 9 
where and Cmatter are given by Eqs. (15.911 and (I5.11|l . and 

Ca{p) '■= ~Ag/i(l + p)^ 6 + (2 + 8P2)p + 2(^2 + Po)p^ 


(A13) 


(AM) 


Eq. (IA13P is also found from Eq. (15.8p when we take the limit of —>■ Ag/3. Hence we conclude that the 
cosmological adiabatic solution discussed in the text corresponds to the present static solution in the “static” limit of 
Ag/3. 


2. Vainshtein screening 

When we assume |A| "C 1, Eq. (IA13I) reduces to 

H^rnls - 2Ag)A + R^pg - K^njpf = 0 , 

which fixes the solution of p by 

„ _ rIPq - K'^R-^Pf 

3(3mgjj — 2Ag) 

From this solution, the ansatz |/t| <C 1 imposes the condition such that 

3|3m2jj - 2Ag| > l^t^pg - K^njpfl. 
In the regime of small p, one of the metric perturbation is given by 




_9 ~ ± 
R ^ 6 


1 + 


m„ 


m„ 


3mgjj 2Ag 1 ® 


RgPg 


3mgjj 2Ag 


^2 ,.2 


Rfpf 




(A15) 


(A16) 


(A17) 


(A18) 


which cannot be restored to the GR result even in the massless limit. 

From Eq. (IA17I) . the ansatz |/t| <C 1 is no longer valid for a high matter density. In the high dense regime, the 
solution of Eq. (IA13|) is given by the equation 




0, 


(A19) 
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from which the typical value of ft is the order of unity. The mass term is negligible compared to the matter term in 
Eqs. (IA8I) - (IA10I) . because KgPg » rrig and KjP/ ^ Wj. Therefore, the metric perturbations in the physical sector 
(and also in the fiducial sector) are restored to the ones in GR. 

However it is not clearly seen whether we can connect the two limiting cases, i.e., the low-dense bigravity region 
with |/1| -C 1 and the high-dense GR region with |/i| ~ 0(1). For the screening, we must find a regular solution 
connecting both regions. For simplicity, we assume that the matter fields are localized in finite regions and only 
consider the outside of the matter distributions. The spatial averages of the matter densities within the radius R and 
Rf are give by 


= PfiR) = 


3Mi 


47tR^(1 + ’ 


where the masses are defined by 


r^o 


Mg = 


AirpgR^dR, M.f = 


rRf,o 


AnpfR^dR , 


(A20) 


(A21) 


respectively. Here Rq and Rf^o are the surface radii of matter distributions. From Eqs. (IA17I) and (IA20I) . the Vainshtein 
radius, beyond which non-linearity of jl becomes important, is estimated as 


^ r 2\GMg-K^gMf\ \ 


1/3 


V - 2A 


(A22) 


We find that the Vainshtein radius depends on the cosmological constant as well as the /-matter field. Outside this 
radius, the smallness condition for p is valid. 

Before introducing the /-matter, we summarize the case without the /-matter as discussed in [^. Since we find 
the similar results in the case with a small mass object of the /-matter, we explain the Vainshtein screening by our 
result. 




FIG. 5: The left figure shows asymptotically homothetic branches (/ —>■ 0) of Eq. (IA13II . We set /32 = —3, bs = 3, = 

rrifjAg = 0, and K^QMf/GMg — 10~® (red solid curve), K^QMf/GMg = 10“^ (blue dashed curve) and K^QMf/GMg = 10 
(green dotted curve). Except for K^QMf /GMg = 10~^, the asymptotically homothetic branches do not extend below the 
Vainshtein radius. The right figure shows the existence of the Vainshtein screening for K^QMf /GMg = 10“®, i.e., GR is 
recovered for R < Rv, where we define = $g|mg=o and 'Fgr = ^ 9 |mg=o. 


In Fig. [5l we present our results in the case of 
K'^QMf = 10~^GMg. The red curves in the figure show 
the value of p (left) and the ratio of the gravitational 
force to the one in GR (right). From these figures, we 
find that GR is recovered within the Vainshtein radius 
Rv^ but the bigravity phase appears in the region of 
R > Rv- This is just the Vainshtein screening. The 


sufficiently small amount of the /-matter does not spoil 
the Vainshtein mechanism. 

We then discuss whether there exists a regular solution 
which shows the Vainshtein screening if we have a large 
amount of the /-matter and/or a cosmological constant. 

First, we consider asymptotically homothetic solutions 
without a cosmological constant. As shown in Fig. [51 if 
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FIG. 6: The same as Fig. [5] for non-homothetic branches. The asymptotically stable non-homothetic branch is characterized 
hy fi ^ 0.70 as R/Rv —>■ oo. The deviation of bigravity from GR diverges as Tg/^'cR oc R^ beyond the Vainshtein radius and 
similar to the $g/$GR. 
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FIG. 7: The roots of Eq. (IA13|) . We set {32 = —3,/I 3 = 
3, rrig = rrij and Ag = lOrrig . An asymptotically stable branch 
are characterized by /i — >■ /loo = 0.28 as R/Rv —>■ 00 . The 
gravitational property is almost same everywhere as discussed 
in Sec.lVH 


the amount of the /-matter held gets large, the /1-curve 
cannot reach below the Vainshtein radius. The deriva¬ 
tive of fi with respect to R diverges around the Vain¬ 
shtein radius, and then the assumption (ITtI) is broken 
there. This hints a sharp transition from the linear to 
nonlinear regimes, and we must not neglect the gradient 
terms of p, such as Rp'^g even if the metric perturbations 
(|l>g|, I'I'gl, |l>/|, and I'F/I) themselves are small. 

There exists another solution which does not approach 
a homothetic solution as R/Rv ^ 1 below the Comp¬ 
ton wavelength of the massive graviton. We call it 
a non-homothetic branch. As shown in Fig. [51 for 
K^QAif/GMg = 10, we find a regular solution which 
approaches a stable non-homothetic branch for R Ry 
(see Appendix lA31 for the stability analysis). In this 
branch, GR is recovered inside the Vainshtein radius, 
but the outside is quite different. The metric pertur¬ 
bations scale as I'bgl, |^g| oc rn^R"^ asymptotically. Note 


that, while the spacetime is given by an inhomogeneous 
spacetime in R , we find the spacetime would ap¬ 

proach an anti-de Sitter spacetime beyond the Compton 
wavelength of the massive graviton by solving the static 
spacetime without the ansatzes (lASp and (lAa [Hi. For 
a smaller amount of the /-matter, however, there is no 
such a regular solution. 

Finally, we consider solutions with a cosmological con¬ 
stant. The behaviour of p is similar to the asymptot¬ 
ically flat case when — 2Ag > 0. However, if 

—2Ag < 0, the behaviour of p changes qualitatively. 
As shown in Fig. [7l when the amount of the g-matter is 
larger than that of the /-matter, the derivative p di¬ 
verges around the Vainshtein radius, indicating a sharp 
transition and the breakdown of our assumption of small 
gradients. On the other hand, we can obtain regular solu¬ 
tions when the amount of the /-matter is larger than that 
of the g-matter. Note that while the poo = 0-branch is 
unstable, the asymptotically poo = 0.280-branch is sta¬ 
ble (see Appendix IA3I for the stability analysis). For 
the stable branch, although GR is recovered within the 
Vainshtein radius, the outside is a stable non-homothetic 
spacetime. As discussed in Sec. El the gravitational be¬ 
haviour is almost the same as in GR when the cosmolog¬ 
ical constant is much larger than the graviton mass. 


3. Stability of a spherically-symmetric static 
spacetime 


We now turn to analyzing the stability of the 
spherically-symmetric static solution obtained in the pre¬ 
vious subsection in the region R ^ in which the 
static solution is also obtained from the cosmological so¬ 
lution ()4.32p in the limit H —>■ 0 or H —>■ Ag/3. The 
matter effects can be ignored in i? Ry- The perturbed 
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metrics are given by 

dsl = + e^^o+'^^odRl + Rldfl^ , (A23) 


as 


dsj = K'^ + R^^'^ 


(A24) 


Tg = T + ST, Rg = R + SR, (A25) 

Tf = T, Rf = {l + fi^)R. (A26) 


where (4>g,'I'g, 4'y) are determined by (IA8|) - (IA11I) 

with /t = /too, a root of 0^^ + Ca = 0. We set the gauge 


J 


We consider the quadratic action in terms of ((/)g,'0g, (5T, 5i?) with AgR^ ^ 1 and rriesR ^ 1 , namely the 

scales smaller than the cosmological horizon and the Compton wavelength of the graviton. We find a constraint 
equation, which is solved when we set 

dlT „ 1 C^TT 


ST =- 


dT 


SR = 


(A27) 


1 + /too l9i? ' 

by introducing the Stiickelberg variable rr. The perturbed metric variables are not dynamical. They are given explicitly 
by 

rniR 


_ ^TT 

i’g = -2^2/ioo + : 


(A28) 


dR 

= 

dR 


(~1 + + (1 + 2/32A00 + 


mlR gjr 

2(l + /too)2(l ^ ^ 


d^n 




1 + 2/ioo + (2/32 - /33)AL 


.527 


(1 + /too)^ 


+ (1 + 2 / 32^00 + P3p‘lc)R'gj,2 


Then we obtain the quadratic action in terms of tt as 


^2 = 


mis 

2tt^ 


J dnj dTdRR^ [z'^'^{dTTrf - Z^^{dRTTf] 


(A29) 

(A30) 

(A31) 


where 




1 


1 + (6/32 — 3) /too + (6/32 ~ 6/32 + 4/33) AL ~ g ~ 10/^3/32 + 2/33 ) AL + ^ AsAL + ^^sAi 


2 r.5 
00 


3 "1/ 


2 1 + Ac 


1 + (6A2 - 2 ) Aoo + (6A2' - 2^2 + 4^3 - 2) AL - 3 (2A2 + (3 - IOA3) P2 - 2^3 + 1 ) A: 
+ 3 (~2A| — 2/32 + A 3 (5^3 + 2)) A^ 


— Ag [1 + ( 4^2 — 2 ) Aoo + (— A2 + 3^3 — 1 ) AL + (As — A2) AL] I 


(A32) 


^RR _ ^ 

“ 2 
3 m 


3 rn^„ 


+ 


2 1 “h jiao 
2 
/ 

2(l + //oo)^ 


1 + (4A2 — 2 ) Aoo + ^ 2 A 2 — 2A2 + AL ~ gAs AL ~ gAs AL 

1 + 4A2A00 + 3 (3A2 + 8A2 + 2A3 — 4) aL + 3 (3A2 ~ A2 + 2A3 — 1) AL 


A„ 


1 + Ao 


+ i (-2A2' + 4 (5A3 - 1 ) A2 - 5 A| + 2A3 - 2) a; 
1+3 (2A2 — 1 ) Aoo ”b f + ^3 ~ 3! i^lc 


(A33) 


r 


The no-ghost and no-gradient instability conditions are Note that for the linear branch Aoo = 0, the coefficients 
given by become 


yRR 

z'^'^ > 0, ^ > 0. 


(A34) 


ZTT = ^ - 2Ag) , (A35) 
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which reproduces Higuchi bound. We also note that for 
Ag, the above expressions reproduce the re¬ 
sults discussed in Sec. II V Cl with w = — 1 by using (14.351) 
and the coordinate transformation. 

For other branches of solutions for jl, since the coef¬ 
ficients are complicated (where /ioo is zero or the roots 
of the sixth-order algebraic equation Cm 2 -I- Ca = 0), we 
cannot obtain the explicit expressions for no-instability 
conditions in terms of /32, /Ja, Wg, m/, and Ag. Hence, we 
check numerically whether the solutions in the nonlinear 
branches are stable or not for some instructive sets of 
coupling constants. 



^2 


FIG. 8: The parameter region of the existence of a stable 
asymptotically non-homothetic branch. We set rUg = 
and Ag = 0. In the colored region, there is a stable branch 
with /Too > — 1 other than jloo = 0. The curves for /Ss = dl 
(solid curve) and Ps = P 2 (dotted line) are also shown for the 
comparison with the stability condition (I4.62|l in the regime 
H 3> uiefi, obtained in Subsec. [IV Cl 


As shown in Fig. |8l the parameters can be classified 
into two types. We show two typical examples in the 
parameter region (14.6211 . 

Case A: We choose 

mg = m/ , ,02 = -3 , /da = 3 , 

in which we can also find a stable branch other than 
Aoo = 0 branch. For instance, the nonlinear stable branch 
is given by 

/too = 0.698872, 

with Ag = 0. For non-zero cosmological constant, the 
stable branch is given by only 

/too = 0.280938, 

with Ag = 10 rUg. 

Case B: We choose 

mg = mf , ^2 = - 1 , ^ , 

in which only /too = 0 is stable with 3mgjj — 2Ag > 0. 
However, although /too = 0 becomes unstable, one stable 
branch appears when — 2Ag < 0. For instance, the 
nonlinear stable branch is given by 

/too = 0.889668, 

with Ag = 10mg. 
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